disp("ProblemC2,C3:");
for n = 1 : 6
    disp("n=");
    disp(n);
    [beta,gamma,alphaP3] = rationalHilbert(n);
    disp("beta=")
    disp(beta);
    disp("gamma=")
    disp(gamma);
    disp("alphaP3=");
    disp(alphaP3);
end

